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. Abstract In brane cosmology, the Big Bang is hypothesized to occur by the anni- 

hilation of the brane-anti-brane pair in a collision, where the branes are three- 
dimensional objects in a higher-dimensional Universe. Spontaneous symmetry 
breaking accompanied by the formation of lower-dimensional topological defects, 
e.g. cosmic strings, is triggered by the so-called 'tachyon condensation', where the 
existence of tachyons is attributable to the instability of the brane-anti-brane sys- 
tem. Here, we discuss the closest analogue of the tachyon condensation in atomic 
Bose-Einstein condensates. We consider annihilation of domain walls, namely 
branes, in strongly segregated two-component condensates, where one compo- 
nent is sandwiched by two domains of the other component. In this system, the 
O ■ process of the brane annihilation can be projected effectively as ferromagnetic 

ordering dynamics onto a two-dimensional space. Based on this correspondence, 
three-dimensional formation of vortices from a domain- wall annihilation is con- 
^ ; sidered to be a kink formation due to spontaneous symmetry breaking in the 

CN , two-dimensional space. We also discuss a mechanism to create a 'vorton' when 

the sandwiched component has a vortex string bridged between the branes. We 
hope that this study motivates experimental researches to realize this exotic phe- 
nomenon of spontaneous symmetry breaking in superfluid systems. 
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1 Introduction 

While a tachyon is a hypothetical superluminal particle in special relativity, a 
tachyon field can exist in quantum field theories, as a consequence of the insta- 
bility of the system. Here, the state considered is initially at a local maximum of 
an effective potential V(T), and then the tachyon 'rolls down' toward a minimum 
of the potential, which causes the exponential amplification of a tachyon field T 
from the initial state T = 0. The rolling-down process is called tachyon condensa- 
tion. In string theory, which is the most promising candidate of the unified theory, 
tachyon condensation occurs in a system containing a Dirichlet(D-)brane and an 
anti-D-brane* 1 . Here, the D-brane is solitonic excitation 2 and the anti-D-brane is 
its anti-object which annihilates the D-brane in collision. The system falls into the 
true vacuum when the annihilation is completed. 

Tachyon condensation of the brane annihilation is applied to brane cosmol- 
ogy^d^ in which phase transitions in the early Universe are supposed to occur 
as a result of a collision of branes. Tachyon condensation causes formations of 
lower-dimensional branes after the collision, which corresponds to nucleation of 
cosmic strings due to the phase transitions accompanied by spontaneous sym- 
metry breaking (SSB) in brane cosmology-^. This phenomenon resembles the 
Kibble-Zurek mechanismiSill, which explains formations of topological defects 
due to conventional SSB in the early Universe 12 . However, SSB due to a brane 
annihilation is different from the conventional SSB in the sense that phase order- 
ing occurs in a lower-dimensional subspace embedded in higher dimensions. The 
understanding of brane annihilation in terms of a phase transition has not been 
yet clarified since how the presence of the extra-dimensions influences the phase 
ordering dynamics is not clear, e.g. the distance or relative velocity between the 
two branes and deformations of the branes into the extra-dimensions. Although 
the Kibble-Zurek mechanism has been tested thoroughly in condensed matter lab- 
oratory experiment s 13 t 14 i 15 i 16 t 17 i 18 i 19 i 20 ^, brane annihilations as sub spatial SSB 
phenomena are quite a novel concept in condensed matter systems. 



Fig. 1 Schematic diagrams of a domain wall annihilation. 

Recently, we showed that vortex formations via the pair annihilation of domain 
walls in two-component Bose-Einstein condensates (BECs) in three-dimensions 
are regarded as kink formations in the effective tachyon field defined in the projected- 
2D space —, as shown schematically in Fig.Q] To the best of our knowledge, this 
system is the first example of subspatial SSB phenomena in condensed matter 
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system. In this system, we can address the nonlinear dynamics of the brane an- 
nihilation, such as defect nucleation and the subsequent effective phase ordering 
dynamics, which is difficult in string theory and brane cosmology. In this paper, 
we briefly review a series of our works on annihilations of domain walls, a 'branes' 
and 'anti-brane', in two-component BEC s 22 i 23 i 24 . 

This paper is organized as follows. In Sec. [2 we introduce the Gross-Pitaevskii 
model and consider a state with a domain wall and an anti-domain wall in strongly 
segregated two-component BEC, where a component is sandwiched between two 
domains of the other component. In Sec. [3j vortex formation from an domain 
wall annihilation is discussed. In Sec. |H the annihilation process is reviewed by 
introducing tachyon field in the projected- 2D conceptually. In Sec. [5] we discuss 
the influence of motion of the sandwiched component along the nucleated vortices. 
Sectionals devoted to a summary and discussion. 



2 Domain wall and anti-domain wall in two-component BECs 

We consider two-component BECs, which consist of condensations of two distin- 
guishable Bose particles. Two-component BECs are well described by the macro- 
scopic wave function as two complex order parameters, *¥\ and *f*>, in the Gross- 
Pitaevskii model at zero temperature 25 . The wave functions obey the action y = 
jdtjd^x^h^dtVi + ^dt^-^-V] with kinetic energy density and 
potential energy density Y\ 

*' = jSll|f|| 4 + jS22| , f!!| 4 + Sl2|f|| 2 | , f!!| 2 -|lll|fl| 2 -^| , f!!| 2 . (2) 

where we used the coupling constant gjk = Inf^a^jm^ (j,k= 1,2) with the 
reduced mass mj^ ^mj 1 -{-m^ 1 and the s-wave scattering length > between 
atoms in the Wj and ^-components. The chemical potential fij > is introduced 
as the Lagrange multiplier for the conservation of the norm Nj = / d 3 x \ *¥j | 2 . Two- 
component BECs undergo phase separation for gyi > ^gngn making domain 
structures with Wj domains (j = 1,2), in which the ^-component is dominant and 
the ^k^j component vanishes. We suppose, for example, the following parameters, 
m\ — M2 = m,g\ \ — g22 = g, = v > an d gn/g = 7=2; this is experimentally 

feasible, e.g. Ref. 26 . We consider the domain structure, where a *f*> domain is 
sandwiched by two *¥\ domains. The length scale of this system is characterized 
by the healing length t, = ft/ \fm\i\ of the ^1 component. 

As a first step of the analysis, we consider the energy density in bulk [the 
potential ^(^i, *f*>) of ©] by neglecting the contribution from the kinetic energy 
(Q]). The amplitude |*£^| is determined to minimize V with W\ fixed. Then the 
density ^2(^1) = l^| 2 * s parameterized by ^1 as n>i — (fi2 — gi2\^i\ 2 )/g for fi2 > 
g^l^il 2 and |^| — for /i 2 < gi2|^i| 2 - By inserting this into y, the potential 
r is reduced to a function W of |*fi|, W^il) = ^(^1, \fn 2 {W{)) = w 4 |^i| 4 + 
^2!^ 1 2 + vt>o- We do not consider the case with v < I/7 since then the potential 
W with W2 < has a local maximum at *¥\ = and *f*> domains cannot exist 
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in bulk. Since the coefficient W4 is negative for /I2 > gi2|^i| 2 but positive for 
M2 < g^l^il 2 * the potential W has a local (global) minimum at l^il = \J\i\jg 
for v > 1 (v < 1). On the other hand, the state with l^l 2 = Pi/g and W\ = 0, 
realized in the *f*> domain, is a local (global) minimum of the bulk potential for 
v < 1 (v > 1). 

Let us consider a pair of flat domain walls, our brane and anti-brane, in an 
uniform system as is illustrated in Fig. [2| (top left). If the thickness of the do- 
main walls is neglected, a stable domain structure is determined from the balance 
between the surface pressures and the interface tension of the walls. The (hydro- 
static) pressure is defined as Pj = gjjnj /2 = fij/2gjj in Wj domain. We impose the 
boundary condition | *Fi (z — >► ±00) | = /g 1 1 , which can be regarded as an external 
pressure P\ — jij/2gjj in this system. A flat domain wall is stabilized for P\ — P2 
with v = 1 . However, for v < 1 , the domain between the two W\ domains can 
be crushed due to the pressure imbalance, P\ > P2, which leads to annihilation of 
the domain walls. This means that, for v < 1, the *f*> domain is destabilized due to 
the existence of domain walls although the domain itself is locally stable far from 
the walls. 

The annihilation should be related to the interaction between the branes. Since 
the 'penetration' of the amplitude l^il (|*£^|) decays exponentially from the do- 
main wall into the *f*> (^1) domain, the short-range interaction between the branes 
is effective only when the inter-brane distance R is comparable to the 'penetra- 
tion depth', called the brane thickness. Here, we define the inter-brane distance 
R as the distance between the two l^i | = |*£*>| planes at z = ±R/2. In our case of 
strong segregated two-component BECs, the brane thickness is of order of healing 
length rsj £ . Here, we focus on drastic annihilation processes caused by the strong 
inter-brane interaction, where the brane and anti-brane are brought close to each 
other at a small R at t = 0. Such a situation can be realized by preparing a stable 
pair of branes with sufficiently large R and reducing rapidly the population of the 
sandwiched component. 



3 Vortex formation from domain wall annihilation 

The conventional process of the pair annihilation is that the brane and anti-brane 
collide straightforwardly to leave only the ^1 component as shown in the ID di- 
agram of Fig. [2] (top). Naively, this process is forbidden because the norm N2 is 
conserved during the process. However, the annihilation process can occur locally 
by pushing away the component. 

In general, the brane annihilation may occur inhomogeneously in various places 
between the branes. Then, the annihilation processes become nontrivial by consid- 
ering the phase difference between the two W\ domains. We introduce the phase 
difference AG with argWi^ = =b°°) = ±AG /2. The annihilation for AG — n is 
just a case, which was realized experimentally by Anderson et al 2 -. In the ex- 
periment—, they created a dark soliton in one component, where the phase of the 
order parameter jumps by n discontinuously across the soliton like a nodal plane, 
and the density dip created by the soliton was filled with the other component. 
Then the filling component was selectively removed, which corresponds to the 
coincident limit R — >> without the *ft component in our model. 
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Fig. 2 Schematic diagram of vortex formations from domain wall annihilation, (top) ID di- 
agrams of domain wall annihilation, (middle) 2D diagrams of vortex formations from domain 
wall annihilation, (bottom) Vortex formation in the 2D diagrams are explained in the configu- 
ration space of W\. In the initial state, the positions P, P\ and P" in the 2D diagram (middle 
left) is located at the origin in the configuration space (bottom left). The direction of motions of 
the points in the configuration space (bottom center) is determined by the directions of super- 
fluid current though each junction (middle center). The circulation along the loop A-B-B"-A"-A 
has no winding number, but the loop for A-B-B'-A'-A has contains a quantized vortex (bottom 
right), whose core is filled with the W2 component (middle right) . 



The inhomogeneous annihilation is initiated by making small junctions, which 
connect locally the two *fi domains by pushing out the component in the region 
—R/2 <z<R/2. Since the phase arg 1 ^ must be continuously connected from a 
domain to the other domain through a junction, the local annihilation causes a 
superfluid current along the z-axis with a current velocity 

hA® „ fiAG-2K ^ 

v z — >0 or v z <0 (3) 

in a m R 

[see Fig. [2] (middle and bottom)]. When the current velocities have the same sign 
between two neighboring junctions (P and P" in Fig. [2]), the annihilation can be 
completed between the junctions. However, if the direction of the velocities are 
different between the two (P and P' in Fig. O, a droplet is left between the junc- 
tion in the 2D diagram and then the phase arg W\ winds once around the droplet, 
making a coreless vortex whose core is filled with the *f*> component. 

The vortex formation strongly depends on how the junctions grow from initial 
fluctuations. For example, in the experiment with small Bose-condensed clouds 
trapped in a spherical harmonic potential—, a vortex ring emerges, affected strongly 
by the geometry of the system. Here, we consider a simplest case of domain wall 
annihilations from initial random fluctuations in an uniform system without exter- 
nal potential. Although the growth rates of junctions with v z > and v z < would 
be different depending on AS, they must be statistically equivalent for AG — n 
because of the symmetry in the configuration space of W\. Especially, if the in- 
stability is sufficiently strong, the fluctuations would grow monotonically and the 
vortex formation would be determined by the distribution of initial fluctuations. 
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Then, a random initial fluctuation develops in a random fashion into the meshed 
structures as shown in Fig. O Since vortices emerge along the boundary between 
the two opposite junctions, the annihilation causes serpentine curves of coreless 
vortices. 




4 Domain wall annihilation viewed in a projected-2D space 

The nucleated vortices stay for a while lying around the z — plane since the z 
component of the velocity, induced by the random distribution of vortices, can- 
cels out on average. In this case, the annihilation process is interpreted as phase 
ordering dynamics of a real scalar field in two dimensions by projecting the three- 
dimensional dynamics onto the x-y plane 22 . When the phenomenon is viewed in 
the projected-2D space, the brane annihilation in the original 3D space may be 
reduced to a ferromagnetic phase transition by a rapid quench. We shall explain 
the correspondence between the 3D and projected- 2D systems qualitatively in the 
following, which is summarized in the table Q] 



3D 


Projected-2D 


Coreless vortices 


Kinks 


Inter-brane distance R 


'Temperature' 


Phase difference A© ^ % 


'External field' 



Table 1 Correspondence between the original 3D system and the projected-2D system. 



Remember that the inhomogeneous annihilation occurs by the growth of junc- 
tions with a superfluid current velocity v z > or v z < 0. In this correspondence, 
we think that the junctions with v z > (v z < 0) is induced effectively by a growth 
of a certain real field T > (T < 0) in the projected- 2D space. Here, the ampli- 
tude and the sign of the field T respectively corresponds to the 'amplitude' and the 
'direction' of the variation of W\'s profile from that in the initial state in the con- 
figuration space as is illustrated schematically in Fig.[2|(see also Fig.0). The state 
of a pair of branes (T = 0) represents the local maximum of the effective potential 
V(T). For the symmetric case with AS — k, the potential should be symmet- 
ric, V(—T)= V(r), and the growth process corresponds to tachyon condensation 
accompanied by spontaneous symmetry breaking, which causes 'magnetization' 
T > or T < 0. 
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The inhomogeneous annihilation sketched in Fig.[3]is projected as an inhomo- 
geneous 'magnetization' onto the projected- 2D space. Since a coreless vortex is 
nucleated between junctions with v z > and v z < 0, a vortex is projected onto the 
2D space as a kink, where the tachyon field T changes its sign across the shadow 
of the vortex line. The tachyon field T is constant, T — ±7&, far from the kink 
core, where the potential V(T) takes the minimum value. Therefore, the domain 
structure of T on a length scale larger than the kink thickness in the projected-2D 
space is independent of how to define the tachyon field T, although the details 
of the internal structure of the kinks depend on the definition of T. It was shown 
numerically that the time development of the domain structure of T obeys actually 
the dynamic scaling law of phase ordering kinetics—. 

The curvature of the potential V(T) at T = is related to the strength of the 
instability, which causes the 'magnetization'. Since the annihilation occurs due 
to the interaction between branes, the instability must become stronger as the 
inter-brane distance R decreases. The curvature of V(T) around T = becomes 
smaller for weaker instability for larger R, and one obtains dV /dT = d 2 V /dT 2 = 
when the inter-brane interaction vanishes completely for R—>oo.Jn terms of the 
Ginzburg-Landau theory for ferromagnetic phase transitions, F2 = d 2 V/dT 2 is a 
function of temperature and becomes negative below a transition temperature. In 
this sense, the inter-brane distance R works as 'temperature'. Since the strength of 
the inter-brane interaction becomes small exponentially with R and the instabil- 
ity should vanish exactly for R —> ©o, the infinity distance may correspond to the 
transition 'temperature'. 

The strength of the instability is related to the line density of vortices nucleated 
immediately after the collision of branes. The line density can be estimated from 
the wave length of the most unstable modes, when the instability is so strong that 
the growth of the mode dominates the vortex formations. According to Ref.— , 
the wave length of the most unstable mode, which has the maximum imaginary 
frequency, increases with R, and thus the line density increases with decrease of R. 
Since the norm N% is an increasing function of R, the nucleated vortices becomes 
thicker as R increases. The /^-dependence of the vortex formation is schematically 
illustrated in Fig.[4j 




Fig. 4 Schematic explanation of the ^-dependence of tachyon condensation, (a) and (b) show 
the cases for smaller and larger R, respectively. The line density of nucleated vortices increase 
with decrease of the initial inter-brane distance R. The nucleated vortices becomes thicker as R 



increases. 
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The potential V(T) becomes asymmetric, V(—T) ^ V(T), for AG ^ n be- 
cause the amplitudes of the superfluid current velocities © induced by the junc- 
tions are different between the two velocities v z > and v z < 0. Similarly to the 
analogy above, the phase difference AG ^ n plays the role of the 'external field' 
in the projected- 2D system. Since the kinetic energy density induced by the su- 
perfluid current is roughly estimated as ~ f yV^, 'magnetization density' T > 
(T < 0) prefers to be amplified energetically for < AG < n (n < AG < 2k). 
Then, the domains of T > (T < 0) tend to percolate across the projected-2D 
system and some domains of T < (T > 0) forms droplets. The droplet cor- 
responds to a vortex ring in the original three-dimensional space. Since a kink 
across the system is realized when both domains (T > and T < 0) percolate, the 
possibility to form vortex lines across the system is maximized for AG = n. The 
A ©-dependence of the vortex formation is schematically illustrated in Fig. [5] 



(a) < A0 < 7T (b) A0 = 7T (c) n < A0 < 2n 




Fig. 5 Schematic explanation of the ^©-dependence of tachyon condensation. (Top) Dynamics 
of tachyon field due to the growth of junctions. Arrows represents the growth of tachyon field 
from the initial state (broken lines) in the configuration space of ^i. The left and right solid 
curves show the mapping of the profiles of the W\ component along junctions with v z < (T < 0) 
and v z > (T > 0), respectively. Here, we neglect the correlation between the two W\ domains 
and assume that the component vanishes between the branes in the initial state with T = 0. 
The domain structure is statistically symmetric between T > and T < domains for (b) A© = 
n, but asymmetric for (a) < AS < % and (c) n < AS < 2k. 



5 Annihilation with bridged vortices between branes 

In the discussion above, we did not take into account explicitly the influence of 
the superfluid current due to the phase gradient of the *f*> component. The influ- 
ence can be important when the population of the *f*> component, localized along 
the coreless vortices, becomes larger. Such a situation is realized in the annihila- 
tion starting from a larger inter-brane distance, where vortices with flat cores are 
nucleated. Especially, when the vortices are bridged between the domain walls, 
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the annihilation dynamics cause an additional nontrivial effect, nucleation of vor- 
ton&2i (see Fig. [6]). 



Fig. 6 Schematic diagram of formation of a vorton from brane annihilation with a bridged vor- 
tex. Solid and broken arrows represent the directions of superfluid currents of the W\ and *f*> 
component, respectively. 

When a vortex ends on the domain wall, the vortex exerts a force on the end 
point joining the brane, causing the wall or the = |*f^| plane to bend. Such a 
configuration is discussed by considering the analogy with string theory 28 , where 
the domain wall and the vortex attached to the wall correspond to a D-brane and 
a fundamental string. If a vortex is bridged between a domain wall and an anti- 
domain wall, the bridged vortex draw both walls [see Fig. [6] (left))]. 

Let us consider a bridged vortex in the ^ component between the two W\ 
domains with a phase difference A & = n. The W\ component has a nodal plane, on 
which we have W\ — 0, and the bridged vortex is perpendicular to the walls passing 
through the nodal plane at the origin O. This configuration is shown schematically 
in the 2D diagram of Fig. Ha) (left). If a junction grows along the bridged vortex 
[Fig- Ha) (center)], two coreless vortices are nucleated below and above the origin. 
Since the component is rotating around the z axis in the initial state, the ¥2 
currents inside the vortex cores flow in opposite directions, positive direction and 
negative direction of the y axis in the 2D diagram of Fig. [2a) (right). 

In 3D diagram of Fig.|7tb), an inhomogeneous growth of junctions can cause a 
vortex loop surrounding the origin. Here, the component inside the vortex core 
has a nontrivial winding along the loop. Such a "twisted" vortex loop is called a 
vorton— In the case of usual vortex loop, in the presence of energy dissipation, 
the loop shrinks easily to smaller loops or rings, and the rings decay into phonons, 
because the energy of a vortex ring is a monotonically increasing function of the 
radius of the ring. However, in the case of a vorton, a "twisted" vortex ring, the 
velocity of the *f*> current along the vortex ring is inversely proportional to the ring 
radius if the winding number of the *f*> component along the ring is fixed. Thus, 
the vorton can be stabilized with a certain radius to reduce the kinetic energy due 
to the ¥2 current. 

The 'twisted' loops can be unstable due to hydrodynamics instability in the 
presence of a relative velocity between the two components. The most likely can- 
didate is the Kelvin-Helmholtz instability, which causes nucleation of vortices 
from the interface between the two phase- separated condensates 29 . The instabil- 
ity may causes vortices in the *f*> component leading to a reduction of the winding 
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Fig. 7 (a) 2D diagrams of brane annihilation with a bridged vortex. A quantized vortices passes 
through the origin along the z-axis. The marks and <g> shows the directions, positive and 
negative directions of the y axis, of the superfluid currents of the *f*> component. Superfluid 
currents of the component are represented with solid arrows, (b) 3D diagrams of the brane 
annihilation viewed from the positive z-axis. Solid and broken arrows represent the directions of 
superfluid currents of the W\ and W2 component, respectively. 

number of the 'twisted' loops. For simplicity, we neglect such an effect and con- 
sider dynamics of the vorton loops in the projected- 2D space by assuming that 
their winding number is conserved during the dynamics. 

Here, we discuss reconnection of the kink loops in the 2D system. A twisted 
loop, along which the *f*> current is winding n times, is called an ^z-loop, where n is 
an integer. An rc-loop with n > (n < 0) is realized, for example, when \n\ bridged 
vortices with a circulation quantum number +1 (—1) are surrounded by the loop. 
Here, we do not consider a kink lying across the system because the *f*> component 
flows along the kink out to the edge of the projected- 2D system, the surface of the 
atomic cloud in the original 3D system, and then the winding number of the kink 
is not well defined. 

(a) rc-loop 72 '-loop (n+n')-loop (b) rc'-loop (n'-n)-\oop 

Fig. 8 Reconnection of twisted loops in the projected-2D space. The integers n, n f , n + n f , and 
n' — n represent the winding number of the W2 current along each loops, (a) A reconnection 
between an w-loop and an n'-loop results in an (« + «') -loop, (b) When an n'-loop surrounds an 
ft-loop, a reconnection between them leads to an (n f — «)-loop. 



A reconnection between twisted loops is classified typically into two cases by 
considering the winding number. The first case is that a reconnection between the 
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two loops, either of which does not surround the other loop [Fig.UJa)]. In this case, 
a reconnection between an rc-loop and an rc'-loop makes an (n + n')-loop. In the 
second case where an rc'-loop surrounds an rc-loop, the reconnection between the 
loops results in an (vl — ft) -loop [Fig. Ob)]. Such nontrivial effects can change the 
statistical properties of the phase ordering dynamics in the projected-2D system. 
For example, the statistic probability of reconnections between loops may depends 
on their winding numbers, which would influence the dynamic scaling law or the 
decay law of kink density in phase ordering kinetics^. 



6 Conclusions and discussion 

Vortex formations via the pair annihilation of domain walls in two-component 
BECs are regarded as kink formations in the effective tachyon field theory of the 
projected- 2D system. The tachyon potential is characterized by the initial inter- 
brane distance and the phase difference between the two bulks separated by the 
brane pair. In the projected-2D system, the inter-brane distance and the phase dif- 
ference play the roles of the temperature and external field in ferromagnetic sys- 
tems, respectively. Vortons can be nucleated via an annihilation of an brane and 
an anti-brane with bridged vortices. 

These phenomena can be monitored directly in experiments of cold atom sys- 
tems with the techniques used in Ref . 27 . The formation of serpentine vortices and 
the effective phase ordering dynamics are feasible by using a trapped condensate 
whose size along the branes is much larger than the characteristic wavelength ~ £ 
of the tachyon mode. The domain structure of the tachyon field can be directly 
visualized by observing the density depletion in the *fi component due to vor- 
tices in the central slice of the atomic cloud—. It is also interesting that a similar 
situations is experimentally realized in superfluid 3 He^2»2I It will be fruitful to re- 
consider the experiments from the viewpoint discussed here. We hope that future 
experiments reveal, for example, how quantum and thermal fluctuations influence 
the defect formations and the relaxation dynamics in the projected- 2D system be- 
yond the mean field theory, which must be useful to establish the SSB phenomena 
occurring in a restricted lower-dimensional subspace. 
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